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Abstract 

A  high  fidelity  computational  campaign  was  initiated  to  predict  impact-induced  detonation.  The 
intent  of  this  numerical  simulation  campaign  is  to  capture  dynamics  and  the  intervening  physical 
processes  leading  to  detonation  of  AIM- 120  Advanced  Medium  Range  Air-to-Air  Missile 
(AMRAAM)  warhead  (WDU-33/B)  resulting  from  fragment  impact.  A  new  and  innovative  first 
principle  calculational  tool  called  “Smoothed  Particle  Hydrodynamics”  (SPH)  is  used.  SPH  is  an 
unstructured  free  Lagrangian  hydrodynamic  tool  developed  earlier  [1,2,3]  for  predicting 
consequences  of  high  strain  rate  reaction  processes  where  materials  experience  extreme  density 
gradients  and  undergo  large  plastic  deformation.  Examples  of  such  impulsive  interactions  can  be 
found  in  penetration  and  fracture/fragmentation  mechanics  problems. 

Well-controlled  field  tests  are  conducted  to  verify  the  detonability  of  the  AIM- 120  class  missile 
warheads  to  high-speed  fragment  impacts.  Detailed  temporally  and  spatially  varying  fragment 
impact  calculations  are  performed  to  illustrate  the  response  of  the  AIM- 120  to  the  kinetic  energy 
of  the  impinging  fragments.  The  correlation  between  calculational  predictions  and  the  field  tests 
data  are  presented. 


Introduction 

A  unique  alternative  to  conventionally  used  structured  Lagrangian  or  Eulerian  analytical  tools  in 
computational  fluid  dynamics  is  Smoothed  Particle  Hydrodynamics  (SPH).  The  SPH  method  is 
an  unstructured,  grid-free,  Lagrangian  scheme  used  to  solve  a  wide  range  of  problems  in 
continuum  mechanics  and  fluid  dynamics.  This  computational  method  simulates  the  motion  of  a 
continuous  medium  by  representing  the  medium  as  an  interacting  finite  system  of  moving 
particles.  This  method  employs  no  grid,  therefore  it  avoids  any  grid  distortion  problems,  which 
occur  in  other  Lagrangian  techniques  when  calculating  large  density  gradients  and  deformations. 
The  principle  of  SPH  was  first  developed  and  introduced  as  a  viable  alternative  to  the  commonly 
used  structured  Lagrangian  method  by  Lucy  [4],  Gingold[5-7],  Monaghan  [8,9],  and  Benz  [10]. 

The  basis  of  SPH  lies  in  the  interpolation  procedure  where  all  thermodynamic  field  variables 
(density,  pressure,  and  internal  energy)  are  evaluated  at  the  discrete  set  of  points  (particles). 
These  points  are  arbitrary  and  randomly  distributed  in  space.  The  interpolation  procedure  uses 
an  interpolating  kernel  function  to  approximate  any  function  of  position  as  an  integral 
transformation.  Using  integration  by  parts,  the  integral  transformation  of  spatial  derivatives  of 
functions  can  be  re-expressed  as  integral  transformations  of  the  function  itself.  These  integral 
interpolants  are  approximated  by  sums  over  a  finite  number  of  particle  positions  or  points.  The 
summation  interpolant  then  replaces  the  governing  partial  differential  equations  of  continuum 
mechanics  with  a  set  of  ordinary  equations. 

Theoretical  development 

The  principles  of  Smoothed  Particle  Hydrodynamics  that  were  introduced  above  will  now  be 
further  described  and  developed.  By  definition,  the  SPH  method  is  based  on  the  motion  of  a 
finite  number  of  interacting  material  particles  distributed  randomly  in  space.  Each  variable  of 
the  continuum  field  in  SPH,  / (r) ,  is  approximated  by  an  integral  interpolant  <  f(r)>  given  by 


1 


(1) 


< 


/ 


(/w)=J/(0 


W 


which  will  be  replaced  by  a  sum  over  the  particle  positions. 


The  function  W  is  the  smoothing  or  interpolation  function  (kernel)  and  has  a  width  measured  by 
h  called  smoothing  length.  For  all  positive  kernels,  i.e.  W>  0,  the  smoothing  function  can  be 
normalized  to  unity. 

,h)  dr  =  l  (2) 

and  it  follows  that  /  (r)  as  /i  — >  0. 

The  mathematical  expression  described  in  (1)  therefore  defines  the  kernel  estimate  (integral 
interpolant)  (/)  of/.  It  is  important  to  note  that  the  smoothing  function  W  resembles  the 
properties  of  Dirac  delta  function  when  h—>  0. 


W  ( 


Suppose  the  function  /  (F)  is  known  only  at  N  discrete  points  (the  particle  positions  rx ,  r2,...), 
i.e.,  they  are  spatially  distributed  according  to  the  number  density  distribution,  rj(r) 


T](r)=^S(r-rj). 

7=1 


Now  if  we  associate  with  the  j-th  particle  a  volume. 


dr ’ 


(3) 

(4) 


we  may  then  approximate  the  integral  interpolant  in  Equation  1  with  the  following  sum  over  the 
particle  positions: 


{f(f)}='Zf(r)w\r-rJ\,h)^-  (5) 

7=1  Pj 


where  =  p (/) . 


The  approximation  methodology  described  above  will  permit  us  to  conveniently  estimate  the 
dynamic  and  thermodynamic  quantities  of  motion  as  described  by  the  Navier-Stokes  Equations 
and  constitutive  models.  Now  in  the  case  where  /(r)  is  replaced  by  the  density  p  (r ) ,  then 
Equation  5  becomes, 

(p(r))  =  ^mjW  (\r-rj\,h)  (6) 

7=1 
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Equation  6  has  a  fundamental  physical  interpretation:  it  states  that  every  material  particle  has  a 
mass  which  is  smoothed  out  in  space  in  accordance  with  the  specific  smoothing  function  W  and 
smoothing  length  h.  Therefore,  density  of  the  medium  can  be  defined  as  a  sum  of  the  density 
distributions  of  all  the  particles  in  the  vicinity  of  the  specific  point. 

1.  The  Kernel,  Smoothing  Function 

There  are  two  interpolation  kernels  or  smoothing  functions  most  commonly  used  in  SPH.  They 
are  Gauss  kernel  and  Cubic  B-Spline  kernel.  For  more  detailed  discussion  see  References  1 1  and 
12. 

2.  SPH  Formulation  of  Hydrodynamics  Equations 

We  use  the  conservation  equations  of  continuum  mechanics  in  the  form: 

Conservation  of  Mass 


dp 

d  t 


pV 


V 


Conservation  of  Momentum 


dV^_=  1  yP 

dt  p 


(7) 


(8) 


Conservation  of  Energy 


E!L=-£^L  vf  va  (9) 

dtp 

Dependent  variables  are  the  scalar  density  (p)  and  the  specific  internal  energy  (E),  the  velocity 

vector  (V) ,  and  the  stress  tensor  o  .  The  independent  variables  are  the  spatial  coordinates 

X  “  and  the  temporal  variable  t.  Here  the  total  time  derivative  ( d/dt )  is  taken  in  the  moving 
Lagrangian  frame.  Summation  over  repeated  Greek  indices  are  implied.  Using  Equation  5,  the 
governing  hydrodynamic  equations  are  approximated  in  the  SPH  framework.  They  become: 

Conservation  of  mass 

-V  j)  •  V*  W^n -rj\ ,  h  )  (10) 
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Conservation  of  momentum 


Conservation  of  Energy 


^r-=2V  £m,fc*-v;)v?wQ;j-r|,fc)  (12) 

d  t  p-  j 

It  should  be  noted  that  an  alternative  form  of  conservation  of  mass  equation  can  be  obtained  from 
Equation  6.  It  reads 


Pt  = 


(13) 


3.  Constitutive  Relations 

The  stress  tensor  appearing  in  the  momentum  and  energy  Equations  1 1  and  12,  respectively,  is 
decomposed  into  an  isotropic  part  which  is  the  pressure  (P)  and  the  traceless  symmetric 
deviatoric  stress  tensor  Sap 


oafi=PSafi-Sa0  (14) 

The  pressure  term  is  calculated  by  means  of  the  equation  of  state  in  the  form  P= p(p,E )  such  as 
the  Mie  Gruneisen  equation  for  solids  or  Gamma  -Law  for  the  gases.  The  two  respective 
equations  are: 

Mie  Gruneisen, 

p(p,E)=(l-0.5rt])PH{p)+TpE 

PH  =ao  n+ba  if  +c0rf  ij>0 

PH=a0J]  ij<0  (15) 

Ideal  Gas  Law, 

P  =  (}-\)pE  0© 

The  subscript  ‘72”  refers  to  the  Hugoniot  curve,  while  7]  =  p  /  p0- 1  represents  the  compression 
and  ris  the  Gruneisen  parameter.  The  constants  a0,  ba,  and  c0  in  Equation  15  are  related  to  the 
coefficients  C  and  S  in  the  shock  and  particle  velocity  relation  Us=C+SUp  through  a  Taylor’s 

series  expansion  of  the  Hugoniot  curve  PH  .  The  ratio  of  heat  capacities  at  constant  pressure  and 
constant  volume  respectively  is  represented  by  }=Cp/Cv. 
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The  anisotropic  part  of  the  stress  tensor  can  be  written  in  a  variety  of  forms  [13,14].  For 
application  in  SPH,  the  deviatoric  stress  field,  assuming  small  displacement,  can  be  written  as 

S“p  =ixlal1  =n(eap  ~Sal>  e”)  (17) 

where  jj,  is  the  shear  modulus  and  £  is  the  strain  rate. 

4.  Artificial  Viscosity 

To  numerically  resolve  sharp  gradients  and  steep  variation  of  the  physical  parameters  at  the 
shock  front,  the  concept  of  artificial  viscosity  is  introduced.  The  artificial  viscosity  function 
smoothes  and  smears  the  shock  front  over  a  few  resolution  lengths.  This  concept  was  first 
introduced  by  Von  Neumenn  and  Richtmyer  [15].  A  special  form  of  the  artificial  viscosity 
function  for  use  in  SPH  was  developed  by  Monaghan  and  Gingold  [7].  It  has  the  following 
form: 


-aCij(Qij  +  B  col  —  —  -  — 

n,  = - H  \  (v;  -Vj)  •  0 ; 

Pij 


n,=o. 


(V-Vj)*(ri-rj)>0 


(18) 


where 


CO;: 


( \ri-rj )2  +eh2 


a~  1.0,/?  *2.0,  £«  0.01 


ft,=l(A  +  ^).C=|(C,  +  C;) 


Here  the  quantity  C  is  the  speed  of  sound. 


Numerical  Simulation 

Detonability  of  AIM- 120  warhead  to  fragment  impact  is  simulated  in  three  dimensions  using  the 
Smoothed  Particle  Hydrodynamics  (SPH)  method  described  in  section  HI.  Broadly  speaking, 
there  are  two  independent  however  interrelated  sequential  phases  to  this  numerical  simulation. 
Each  phase  is  simulated  in  a  continuous  numerical  process  which  applies  fundamental  physics  to 
assess  the  complex  interactions  between  the  two  phases.  The  first  phase  of  simulation  consists  of 
detonation,  breakup  and  fragmentation  of  the  AIM-120  donor  warhead.  The  importance  of 
solving  the  donor  explosion  from  a  first  principles  foundation  is  apparent.  One  must  accurately 
characterize  fragment  shape  and  distribution,  incident  impact  angle,  and  impact  speed  before 
predicting  failure  modes  (i.e.,  detonability)  of  the  warhead  acceptor.  In  this  simulation  some 
effects  such  as  fragment  distortion,  inter-fragment  distance,  rotation,  or  skewness  are  considered 
second  order  effects  and  are  ignored. 
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Next,  propagation  of  the  radially  expanding  fragments  is  modeled  and  tracked  both  spatially  and 
temporally  before  they  impact  a  second  AIM- 120  warhead,  referred  to  as  the  acceptor  warhead. 
The  acceptor  warhead  is  positioned  at  a  specified  configuration  and  at  a  prescribed  standoff 
distance  with  respect  to  the  donor  warhead.  The  second  phase  of  the  simulation  consists  of 
determining  the  response  of  the  warhead  acceptor  to  the  kinetic  energy  of  the  impacting  donor 
fragments.  The  complex  phenomena  involved  in  the  impact  process,  including  possible 
perforation  of  the  warhead  casing  and  the  ensuing  shock  loading  of  the  PBX  (AF)-108  explosive, 
is  simulated  at  this  stage.  Finally,  the  induced  pressure  level  in  the  acceptor  warhead  is  assessed 
to  determine  whether  it  exceeds  the  threshold  value  required  for  sustaining  a  detonation.  The 
threshold  pressure  value  necessary  for  initiating  a  self-sustaining  detonation  front  in  PBX  (AF)- 
108  is  estimated  to  be  20  kbar. 

1.  AMRAAM  Warhead  Geometry  Setupf 

The  physical  geometry  of  the  donor  and  acceptor  warheads  is  constructed  (generated)  in  three 
dimensions  using  the  SPH  code  MAGI.  The  AISI-  4140  medium  carbon-chromium  steel  used  in 
the  warhead  casing  is  modeled  by  the  Mie-Gruneisen  Equation  of  State  and  the  Elastic-Perfectly- 
Plastic  constitutive  strength  model.  Special  coding  was  developed  to  construct  interior  and 
exterior  scoring  of  the  warhead  casing.  To  capture  dynamics  of  the  case  rupture  along  the  score 
lines  and  fragmentation,  the  Johnson  and  Cook  fracture  model  [16]  was  used.  The  Johnson- 
Cook  fracture  model  uses  a  failure  criteria  based  on  equivalent  plastic  strain,  taking  into  account 
the  pressure,  temperature,  and  strain  rate  along  the  loading  path  for  each  material  element.  The 
warhead  model  with  the  explosive  removed  is  shown  in  Figure  1.  The  hanger  assembly  shown 
in  this  figure  was  placed  in  the  model  in  anticipation  of  assessing  the  resulting  fragment 
distribution  and  its  effect  on  an  acceptor  warhead.  For  the  simulations  reported  here,  the  hanger 
was  oriented  away  from  the  acceptor  warhead.  Further  specifics  of  the  warhead  geometry  and 
problem  setup  are  given  in  Figure  2.  This  3-D  model  was  constructed  from  warhead  engineering 
drawings  and  used  909,578  particles  resulting  in  a  resolution  of  1.7  mm  (particle  diameter)  in  the 
acceptor  explosive  material  and  the  donor  and  acceptor  casing.  This  resolution  was  verified  in  a 
2-D  convergence  analysis  described  later.  The  3-D  simulation  used  two  reflective  planes  of 
symmetry  along  the  donor  and  acceptor  warhead.  The  simulations  were  executed  on  a  shared 
memory  parallel  machine  and  required  about  1000  cpu  hours  among  16  processors. 


f  calculations  are  performed  by  Mr.  David  Medina  at  AFRL/DEPA. 
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Figure  2.  Problem  setup  parameters. 
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2.  Explosive  Reaction  Simulation 

A  total  of  6783  grams  equivalent  PBX  (AF)-108  explosive  was  used  in  the  simulation  of  the 
AMRAAM  (AIM-120)  warhead  (WDU-33/B).  It  is  important  to  note  that  the  6783  grams  PBX 
used  in  the  simulation  consisted  of  two  parts;  6652  grams  main  charge  plus  131  grams 
equivalent  PBX  to  account  for  the  1 10  grams  of  CH-6  explosive  used  in  the  warhead  booster. 
This  energy  equivalent  conversion  is  based  on  a  1.189  conversion  factor  obtained  from  the  ratio 
of  mass  density  of  RDX  used  in  CH-6  to  that  of  PBX. 

The  unreacted  solid  explosive  used  in  these  simulations  was  represented  by  the  Mie-Gruneisen 
equation  of  state.  The  gas  phase  of  the  reacted  explosive  was  modeled  with  the  JWL  equation 
of  state  (EOS).  The  JWL  EOS  describes  the  pressure,  volume,  and  energy  relationships  in  the 
adiabatic  expansion  region  of  the  explosive  products,  given  by 


P  =  A 

1  a 

e-Ry  +  B 

RV 

l  a  e~w  | 

R,V  V 


(19) 


For  isentropic  processes,  the  pressure  (Ps)  is, 


Ps  =  AeRy  +  BeRy  +  CV“(a>+1) 


(20) 


where  P  is  the  pressure  of  the  detonation  products,  A,  B,  and  C  are  linear  coefficients,  Ri,  R2, 
and  (0  are  nonlinear  coefficients,  V  is  the  ratio  of  the  detonation  products  volume  (v)  to  the  initial 
explosive  volume  (v0),  and  E  is  the  detonation  energy  per  unit  volume. 

The  physical  constants  specific  to  PBX  (AF)-108  were  obtained  from  References  18,19,  and  20. 
These  values  are  summarized  in  Tables  1  and  2.  Other  relevant  input  data  to  the  MAGI 
constants  file  is  given  in  Table  3.  It  should  be  noted  that  the  detonation  threshold  (item  78  in 
constants  file)  was  obtained  from  Wagenhals  [21]. 


TABLE  1.  PBX  (AF)-108  EQUATION  OF  STATE  PARAMETERS 


Density 

(g/cm3) 

Dv 

(cm/psec) 

E19 

(kJ/g) 

Pcj 

(kbar) 

Gamma 

Y 

AH(je|  1 

(Mbar-cc/g) 

1.57 

0.791 

1.095 

260 

2.784 

0.0453 

TABLE  2.  JWL  PARAMETERS  FOR  PBX  (AF)-108. 


vch 

Rl 

R2 

CO 

A(Mbar) 

B(Mbar) 

C(Mbar) 

0.7358 

4.12 

1.031 

0.3 

4.725 

0.03230 

0.01179 
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The  Forest  Fire  Reactive  bum  model  is  used  to  calculate  the  incipient  detonation  in  the  acceptor 
warhead.  In  this  reactive  bum  model,  the  results  of  simple  small-scale  wedge  tests  are  applied  to 
calibrate  the  model,  allowing  one  to  extrapolate  to  complex  problems  involving  shock-initiated 
reaction  of  heterogeneous  solid  explosives.  The  derivation  of  this  method  is  based  on  the 
assumption  of  a  global  reaction  rate  model  for  the  decomposition  of  the  explosive.  The  model 
has  been  successfully  applied  to  a  variety  of  very  practical  problems.  For  a  detailed  derivation, 
the  reader  is  referred  to  Mader  [22]. 

TABLE  3.  RELEVANT  VALUES  IN  MAGI  CONSTANT  FILE. 


afx 

!  pbx  {af)-108 

□ 

1 

1 . 6500E-01 

co  -  ambient  sound  speed  (cm/us)  Navy  Expl  HB,  Pg  4-24 

2 

1.9600E+00 

s  -  slope  of  us  -  up  line  Navy  Expl  HB,  Pg  4-28 

11 

1.5000E+00 

Gamma  -  Gruneisen  constant 

12 

2 . 5795E-01 

Cv  -  specific  heat  (cal/(g-K)  Navy  Expl  HB,  Pg  4-33 

13 

6 . 3776E-01 

Vo  -  ambient  specific  volume  (cm* *3}  Navy  Expl  HB,  Pg  1-5 

14 

7 . 0000E-05 

alpha  -  coefficient  of  linear  expansion,  Navy  Expl  HB,  Pg  4-35 

15 

2 . 5000E-01 

mu  -  shear  modulus  (mb) 

16 

3 . OOOOE-03 

Yo  -  yield  strength  (mb) 

17 

3 . OOOOE+02 

To  -  ambient  temp  (K) 

18 

2 . 5000E+03 

Tmelt  -  melt  temp  (K) 

19 

-7 . 03  00E-06 

Pmin  -  spall  strength  (mb)  Navy  Expl  HB,  Pg  5-11 

20 

2 . 0000E-Q2 

Multiplier  for  density  clip 

68 

2 . 2213E+02 

mw  -  molecular  weight  for  RDX  (72%  AFX)  NEH  Pg  2-8 

69 

2 . 7848E+00 

gamma  -  specific  heat  ratio 

70 

4 . 7250E+00 

jwl  a  -  Eglin  Energetic  Materials  Branch 

71 

3 . 23  00E-02 

jwl  b  -  Eglin  Energetic  Materials  Branch 

72 

1.1790E-02 

jwl  c  -  Eglin  Energetic  Materials  Branch 

73 

4 . 1240E+00 

jwl  rl  -  Eglin  Energetic  Materials  Branch 

74 

1 . 03 10E+00 

jwl  r2  -  Eglin  Energetic  Materials  Branch 

75 

3 . OOOOE-Ol 

jwl  w  -  Eglin  Energetic  Materials  Branch 

76 

7.9170E-01 

detonation  velocity  -  cm/us 

77 

4 . 53  84E-02 

detonation  energy  -  mb-cc/gm  Eglin  Energitic  Mat 

78 

1 . 8000E-02 

detonation  threshold  for  burn  2  (currently  mb) 

79 

2 . 6000E-01 

Chapman- Jouguet  detonation  pressure,  From  Eglin 

Simulation  Results 

The  susceptibility  of  the  AMRAAM  (AIM- 120)  warhead  to  fragment  impacts  resulting  in 
possible  sympathetic  detonation  was  simulated  for  two  distinct  standoff  distances.  The  donor 
and  the  acceptor  warheads  were  located  at  5  and  33.5  inches  standoff  distances  from  each  other 
in  a  perfect  center-to-center  alignment  configuration.  Following  propagation  of  the  donor 
warhead  fragmentation,  it  becomes  apparent  that  the  likelihood  of  perfect  normal  impact  on  the 
acceptor  warhead  decreases  with  increasing  cylindrical  expansion  of  the  fragments  (Figure  18). 
The  current  3-D  simulations  take  this  expansion  into  account  over  the  calculational  process. 

1.  Diagnostic  Calculations 

A  set  of  preliminary  calculations  were  conducted  to  evaluate  the  properties  of  the  various 
physical  modules.  Additionally,  attempts  were  made  to  assess  the  convergence  of  the  results  as  a 
function  of  the  resolution  particle  size  and  setup.  A  2-Dimensional  Cartesian  coordinates  model, 
shown  in  Figure  3,  was  used  to  model  strip  fragment  generation  and  its  propagation. 
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The  smallest  size  fragment  allowed  in  this  analysis  was  a  fragment  formed  by  the  scoring 
pattern,  i.e.,  14.8  x  11.7  mm2. 

To  ensure  that  the  scored  casing  was  rupturing  properly  and  the  brittle  fracture  module  behaved 
realistically,  we  compared  our  calculated  rupture  radius  with  the  predicted  rupture  radius,  for  an 
unscored  casing,  obtained  from  the  Gurney  Equation,  reference  23.  As  expected,  the  calculated 
scored  casing  rupture  occurs  at  a  reasonable  time  before  the  unscored  rupture  (Figure  4). 

MATERIAL  PLOT 
TIME  =  0.00  MICRO-SEC 


F  i  .  ""’i . r  i  i  . . t . t . — i - * 

15.00  30.53 


Figure  3.  Two  Dimensional  Model. 
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Original  radius.  Scored  rupture  Rupture  radius  for  unscored  AISI-4140  casing 

Ro  =  8.89  cm  occurring  about  R  =  16  cm.  R/R0  =  1.8 

12  cm 

Figure  4.  Rupture  characteristics  of  the  WDU-33/B  Warhead. 

Furthermore,  the  acceleration  and  the  deceleration  character  of  the  donor  fragments  were 
assessed  separately  using  a  2-D  analysis.  In  this  analysis,  it  was  assumed  that  all  the 
aerodynamic  forces  including  drag  and  friction  forces  were  negligible.  The  2-Dimensional 
calculation  showed  that  fragments  reach  a  steady-state  velocity  at  about  90  microseconds 
(Figure  5). 

2-D  Fragment  Expansion  Velocity 


It  is  worth  noting  that  a  threshold  impact  velocity  for  detonation  was  established.  The  lower 
bound  of  impact  velocity,  which  did  not  detonate,  was  3608  ft/s  while  the  upper  bound,  which 
resulted  in  detonation,  was  3973  ft/s  (Figure  6). 
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Finally,  the  threshold  impact  velocity  for  a  single  WDU-33/B  fragment  impacting  the  same 
warhead  was  calculated  using  the  Jacobs-Roslund  empirical  formula 

Vc  =(A/  dl/2)(l+ B)(\+CT/  d)  (21) 

where  Vc  =  critical  impact  velocity  for  target  detonation 

d  =  fragment  critical  dimension,  e.g.  diameter  (mm) 

T  =  target  cover  thickness  (mm) 

A  =  explosive  sensitivity  coefficient  (mm3/2/psec) 

B  =  fragment  shape  coefficient  (dimensionless) 

B  =  0  for  flat-end  fragments 
C  =  cover  plate  protection  coefficient  (dimensionless) 

The  values  used  for  PBX  (AF)-108  for  A,  B,  and  C  are  given  in  the  Navy  Explosives  Handbook 
[24],  From  this  formula,  the  threshold  velocity  for  a  WDU-33/B  fragment  impacting  the  same 
warhead  is  6876  ft/s. 

There  is  a  tight  connection  between  the  physical  size  of  the  impacting  fragments,  failure,  or 
critical  diameter  of  the  explosive  and  the  computational  resolution.  The  critical  explosive 
diameter  is  the  minimum  diameter  at  which  a  cylindrical  charge  of  HE  sustains  a  high  order 
detonation.  Although  there  exists  adequate  information  on  critical  diameter  of  RDX  based 
explosive,  none  was  found  for  PBX  (AF)-108  in  References  24  or  25. 

As  shown  in  Figure  7,  the  individual  fragment  size  in  the  WDU-33/B  Warhead  follows  the 
scoring  dimensions  and  is  most  likely  well  above  the  critical  diameter  of  the  PBX  (AF)-108 
explosive.  Therefore,  to  ensure  high  resolution,  a  computational  (SPH)  particle  size  of  1.67  mm 
was  selected  for  this  class  of  impact  induced  detonation  problem. 


12 


Figure  6.  Two-dimensional  calculation  of  the  WDU-33/B  casing  being  impacted 
with  a  warhead  fragment  at  (a)  3608  ft/s  and  (b)  3973  ft/s. 


40.1762  mm 


f _ 


11.684  mm 

T- 


14.7502  mm 


Figure  7.  Dimensions  of  a  single,  undeformed  fragment  in  the  WDU-33/B  warhead. 


2.  Five  Inch  Standoff  Results 

A  standoff  distance  of  five  inches  was  selected  as  a  lower  bound  for  simulating  WDU-33/B 
warheads  in  alignment  with  each  other.  The  problem  setup  is  illustrated  in  Figure  2.  The  donor 
warhead  was  initiated  in  the  aft  (top)  center  position  where  the  fuse  would  be  located.  The 
position  where  detonation  is  initiated  is  important  because  it  controls  the  fragment  trajectories  as 
shown  below.  For  the  five-inch  standoff,  the  condition  of  the  casing  just  prior  to  contact  with  the 
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acceptor  (100  microseconds)  is  shown  in  Figure  8.  The  hanger  assembly  shown  in  this  figure 
was  not  oriented  to  impact  the  acceptor  since  this  would  be  the  configuration  for  a  wing  mounted 
missile.  As  expected,  the  rings  of  free  fragments  around  the  center  have  the  highest  radial 
expansion  speed  outward.  The  top  and  bottom  closure  plates,  along  with  the  top  and  bottom 
rings  of  fragments,  are  maintained  as  a  contiguous  mass  for  a  longer  time  since  the  top  of  these 
fragment  rings  did  not  have  a  score  line. 


Figure  8.  Early  time  explosive  expansion  of  the  donor  warhead. 


The  evolution  of  the  explosive  expansion  can  be  seen  in  the  2-D  slice  of  the  3-D  simulation. 
This  is  shown  in  Figure  9.  At  the  time  impact  occurs  (100  microseconds),  rupture  of  the  donor 
case  is  starting,  so  impact  on  the  acceptor  is  closer  to  that  of  case  slap.  The  detonation  point  of 
the  donor  warhead  along  the  top  center  has  resulted  in  a  pressure  profile  that  induces  the  center 
fragments  to  move  slightly  downwards.  This  motion  becomes  more  obvious  for  the  33.5  inch 
standoff. 
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MATERIAL  PLOT  MATERIAL  PLOT 

TIME  =  100.00  MICRO-SEC  TIME  =  110.00  MICRO-SEC 


0.00  32.00  0.00  32.00 


Figure  9.  Two-dimensional  slice  showing  evolution  of  the  donor  warhead  expansion. 

Detonation  of  the  acceptor  is  underway  by  about  110  microseconds  (Figure  10).  The  peak 
pressure  is  300  kbar  and  is  that  of  the  detonation  front  in  the  explosive. 


Figure  1 1  is  a  view  from  the  position  of  the  donor  warhead  and  shows  the  incoming  (outgoing 
from  the  view  location)  fragments  as  they  approach  the  acceptor  warhead.  The  view  is  from 
inside  the  donor  warhead  looking  towards  the  acceptor.  By  this  time,  all  chemical  energy  (30.8 
M Joules)  is  liberated.  In  Figure  12,  the  deformation  of  the  acceptor  outer  casing  is  shown  at  120 
microseconds.  At  this  time,  detonation  is  underway  but  has  not  had  time  to  cause  the  casing  to 
visibly  expand.  The  expansion  of  the  warhead  casing  due  to  detonation  is  clearly  evident  by  140 
microseconds  (Figure  1 3). 

PRESSURE  (KB) 

TIME  =  110.00  MICRO-SEC 


3.7GE+02 

3.G8E+02 

2.46E+02 

1  . 84E+C2 

1  . 23E+02 

6 , 1 GE+0 1 

-6.73E-01 

-6  ►  24E+0 1 
-1  >  24E+02 
-1  .SftE+02 

-2.4SE+02 

0.00  32.00 


Figure  10.  Pressure  state  of  acceptor  warhead  at  incipient  detonation. 
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Figure  11.  Case  expansion  shortly  after  impact  (110  microseconds). 


Figure  12.  Acceptor  casing  damage  for  5  inch  standoff  distance  at  120  microseconds. 
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Figure  13.  Acceptor  casing  showing  explosive  expansion  for  5  inch  standoff  at  140  psec. 

3.  33.5  Inch  Standoff  Calculation 

A  standoff  distance  of  33.5  inches  was  selected  for  this  simulation  with  the  intention  of 
replicating  the  standoff  distance  between  AIM- 120  missiles  in  the  second  and  third  wing 
mounting  positions  (the  first  position  being  the  wing  tip  position)  of  an  F- 1 6  aircraft.  Hand 
measurements  were  taken  from  a  fully  loaded  F-16  (Figure  14).  The  relative  offset  of  the  two 
missiles  along  the  missile  axis  (measured  to  be  about  27  inches)  was  not  included.  As  noted 
previously,  and  in  coordination  with  ongoing  testing  program,  the  warheads  in  these  simulations 
were  in  direct  alignment  across  each  standoff  distance. _ 


Figure  14.  F-16  wing  mounted  AIM- 120s. 
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The  problem  setup  was  similar  to  that  for  the  previous  simulation  except  for  the  standoff  distance 
(Figure  2). 

Since  the  hanger  assembly  was  in  the  vertical  upward  position  for  both  wing  mounted  missiles,  it 
would  not  be  ejected  in  the  direction  of  the  adjacent  missile  in  a  detonation.  For  this  reason,  it 
was  not  included  in  this  simulation.  The  larger  standoff  distance  in  this  simulation  offered  more 
opportunity  for  the  fragment  distribution  to  be  affected  by  the  initial  conditions.  For  example, 
the  slight  downward  component  of  velocity  for  the  center  fragments  caused  the  lower  ring  of 
fragments  to  narrowly  miss  the  acceptor  (Figure  15).  This  2-D  slice  was  taken  at  3.79  cm  from 
the  centerline.  In  this  figure,  the  fourth  fragment  from  the  bottom  disappears  by  550  ps  because 
it  moves  out  of  the  2-D  slice  plane.  Also  resulting  from  the  long  standoff  distance  is  the  wide 
separation  between  the  radially  expanding  fragments.  By  the  time  they  impact,  they  have 
separated  enough  for  the  acceptor  warhead  to  fit  through,  had  the  angular  position  of  the  donor 
been  adjusted  by  a  few  degrees  prior  to  detonation.  For  the  initial  orientation  of  the  donor  in  this 
simulation,  the  resulting  radial  trajectory  of  the  fragments  and  the  final  impact  position  is  shown 
in  Figures  16  and  17.  The  radial  position  of  the  fragment  first  to  impact,  and  the  one  adjacent  to 
it,  is  shown  in  the  2-D  slice  view  shown  in  Figure  17.  As  seen  in  this  figure,  the  initial  angular 
position  of  the  score  marks  in  the  donor  warhead  is  critical  to  the  location  of  impact  on  the 
acceptor  warhead.  A  few  degrees  rotation  in  the  counter-clockwise  direction  of  the  donor 
warhead  would  have  produced  a  fragment  distribution  that  missed  the  acceptor  warhead. 

From  the  pressure  contours  shown  in  Figure  1 8,  it  can  be  seen  that  detonation  is  initiated 
primarily  by  the  shock  induced  by  the  second  fragment  from  the  bottom. 

The  three  dimensional  view  of  the  donor  fragments  at  impact  is  shown  in  Figures  19  and  20.  In 
Figure  20,  explosive  expansion  is  apparent  in  the  casing  by  570  microseconds.  The  plastic 
deformation  in  the  acceptor  casing,  due  primarily  from  fragment  impact,  is  shown  in  Figure  21. 

=  520.00  MICE  530.00  MICE!  540.00  MICRO  550.00  MICR‘ 


Figure  15.  Progression  of  impacting  fragments  on  the  acceptor  warhead 

with  a  33.5  inch  standoff. 
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530  micro-sec 


-20.00  20.00 


550  micro-sec 


Figure  16.  Radial  distribution  of  fragments  at  impact  for  the  33.5  inch  standoff. 
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Figure  17.  Fragment  evolution  for  the  33.5  inch  standoff  simulation. 
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Figure  20.  Acceptor  warhead  at  570  psec  showing  explosive  expansion. 


Figure  21 .  Damage  to  acceptor  casing  at  the  start  of  induced  detonation  (540  microseconds). 
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Experimental  Effort 


Well-controlled  and  highly  instrumented  full-scale  field  experiments  will  be  conducted  to  assess 
susceptibility  of  the  AIM- 120  warhead  to  fragment  impacts.  There  are  two  main  objectives: 

1 .  To  experimentally  investigate  probability  of  fragment  impact  induced  detonation  of  the 
AIM- 120  warhead  (WDU-33/B)  and  occurrence  of  Maximum  Credible  Event  (MCE) 
as  a  function  of  standoff  distance,  and 

2.  To  provide  high  fidelity  fragmentation  and  impact  data  essential  for  calibration  and 
validation  of  the  analytical  tools  used  for  characterizing  impact  induced  detonation 
phenomena. 

A  myriad  of  diagnostics  including  Time  Of  Arrival  (TO A)  pins,  pressure  gages,  and  witness 
plates  are  placed  at  prescribed  ranges  around  the  test  bed  to  measure  emanating  shock  wave 
velocity  and  pressure  as  well  as  monitoring  fragments  propagation.  Other  transducer  and 
monitoring  systems  used  are  Photodiodes  and  high-speed  cameras.  Photodiodes  are  placed  on 
the  outer  surface  of  the  donor  warhead  to  detect  first  light  break-out  of  the  warhead  indicating 
warhead  case  rupture.  The  high-speed  cameras  will  be  used  to  record  formation  of  the  fireball. 

1.  Test  Plan 

Under  the  auspices  of  the  Defense  Threat  Assessment  Agency  Albuquerque  Office  (DTRA  AO), 
Applied  Research  Associates,  Inc.  has  developed  a  comprehensive  test  plan  to  conduct  six  AIM- 
1 20  warhead  detonation  tests  [27] .  Figures  22  through  26  illustrate  specifics  of  these  test  series 
including  test  bed  layouts,  warhead  donor/acceptor  placement  and  the  instrumentation 
configurations. 

It  should  be  noted  that  no  test  was  done  while  this  paper  was  being  prepared  and  therefore,  no 
correlation  with  the  calculational  prediction  was  performed. 
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Conclusions 


A  novel  first  principle  particle  mechanics  code  ( SPH )  was  used  to  simulate  the  dynamics  of 
fragment  impact  induced  detonation  of  the  WDU-33/B  warhead.  Two  distinct  standoff  distances 
of  5  and  33.5  inches  were  used  in  this  simulation.  An  ideal  configuration,  where  the  acceptor 
and  the  donor  warheads  were  in  perfect  alignment,  was  assumed  in  these  calculations. 

For  these  special  ideal  cases  where  aerodynamic  effects  were  not  included,  simulation  predicted 
acceptor  warhead  detonating  at  both  standoff  distances.  It  was  shown  that  the  likelihood  of  a 
perfectly  normal  fragment  impact  on  the  acceptor  warhead  decreases  with  the  increasing  standoff 
distance  which  consequently  diminishes  the  likelihood  of  sympathetic  detonation.  Furthermore, 
fragment  distribution  pattern  indicates  that  the  introduction  of  a  relative  axial  offset  between  the 
donor  and  acceptor  warheads  will  greatly  influence  deposition  of  kinetic  energy  in  the  acceptor 
warhead  and  hence  the  possibility  of  sympathetic  detonation.  Exact  details  of  these  geometrical 
and  physical  phenomena  will  be  fully  explored  once  the  experimental  results  become  available. 

Finally,  as  has  been  demonstrated,  dynamics  of  the  warheads  interactions  under  high-strain  rate¬ 
loading  due  to  the  high  velocity  fragment  impact  and  all  the  intervening  processes  including 
warhead  rupturing,  fragmentation,  and  fragment  distribution  can  be  simulated  numerically. 
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